The Tunka-133 EAS Cherenkov light array: status of 2011 
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Abstract 

| A new EAS Cherenkov light array, Tunka-133, with ~1 km 2 geometrical area has been installed at the Tunka Valley (50 km from 
["t"] Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 
jTj '10 16 - 10 18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration. The method 
of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS 
Cherenkov light allows to estimate the depth of the EAS maximum X max . 

The results on the all particles energy spectrum and the mean depth of the EAS 
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lx vs. primary energy derived from 



O the data of two winter seasons (2009 - 2011), are presented. Preliminary results of joint operation of the Cherenkov array with 
antennas for detection of EAS radio signals are shown. Plans for future upgrades - deployment of remote clusters, radioantennas 
^2 and a scintillator detector network and a prototype of the HiSCORE gamma-telescope - are discussed. 

Keywords: EAS Cherenkov light array, cosmic rays, energy spectrum and mass composition. 
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1. Introduction 

The study of primary energy spectrum and mass composi- 
tion in the energy range 10 15 - 10 18 eV is of crucial importance 
for the understanding of the cosmic rays origin and propagation 
in the Galaxy. 

To measure the primary cosmic ray energy spectrum and 
mass composition in the mentioned energy range, the new array 
Tunka-133 (HI], J2l), with nearly 1 km 2 geometrical area has 
been deployed in the Tunka Valley, Siberia. It records EAS 
Cherenkov light using the atmosphere of the Earth as a huge 
calorimeter and has much better energy resolution (~ 15%) than 
EAS arrays detecting only charged particles. 

2. Tunka-133 

The Tunka-133 array consists of 133 wide-angle optical de- 
tectors based on the PMT EMI 9350 with a hemispherical pho- 
tocathode of 20 cm diameter. The detectors are grouped into 
19 clusters, each cluster with seven detectors - six hexagonally 
arranged detectors and one in the center. The distance between 
the detectors inside the cluster is 85 m. 

The Cherenkov light pulse is sent via 95 m coaxial cable 
RG58 to the center of a cluster and digitized. The dynamic 



range of the amplitude measurement is about 3 • 10 4 . This is 
achieved by means of two channels for each detector extracting 
the signals from the anode and from an intermediate dynode of 
the PMT with different additional amplification factors. 

The cluster electronics includes the cluster controller, 4 four- 
channel FADC boards, an adapter unit for connection with op- 
tical dtectors and a special temperature controller. The 12 bit 
and 200 MHz sampling FADC board is based on AD9430 fast 
ADCs and FPGA XILINX Spartan XC3S300 microchips. The 
cluster controller consists of an optical transceiver, a synchoniza- 
tion module, a local time clock and a trigger module. The op- 
tical transceiver operating at 1000 MHz is responsible for data 
transmission and formation of 100 MHz synchronization sig- 
nals for the cluster clocks. The cluster trigger (the local trigger) 
is formed by the coincidence of at least three pulses from the 
optical detectors exceeding the threshold within a time window 
of 0.5jus. The time mark of the local trigger is fixed by the clus- 
ter clock. The accuracy of the time synchronization between 
different clusters is about 10 ns. Each cluster electonics is con- 
nected with the DAQ center by a multi-wire cable containing 
four copper wires and four optical fibers. 

The central DAQ station consists of 4 DAQ boards synchro- 
nized by a single 100 MHz generator. The boards are connected 
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Figure 1 : : Example of a pulse from one Tunka- 133 detector, recorded by high 
and low gain channels 

with the master PC by 100 MHz Ethernet lines. 

3. Data processing and reconstruction of the EAS parame- 
ters 

3.1. Data processing 

The primary data record for each Cherenkov light detector 
contains 1024 amplitude values with step of 5 ns (FigQ]). Thus 
the each pulse waveform is recorded together with the preced- 
ing noise as a trace of 5 fis. To derive the three main parameters 
of the pulse: front delay at a level 0.25 of the maximum ampli- 
tude fj, pulse area Q; and full width at half-maximum (FWHM) 
t,, each pulse is fitted with a special self-designed smoothing 
curve f^. 

The reconstruction of the EAS core position is performed 
with two methods - by fitting the measured charges Q, with the 
lateral distribution function (LDF ([4]) and by fitting the mea- 
sured pulse widths r, by the width-distance function (WDF)(|5]). 

3.2. Energy reconstruction 

As a measure of energy we used the Cherenkov light flux 
density at a core distance of 175 m - 2(175). Connection be- 
tween the EAS energy Eo and Q175 may be expressed by the 
following formula: 

Eo = C ■ Q\ 1S (1) 

It was found from CORSIKA simulation, that for the en- 
ergy range of 10 16 - 10 17 eV and zenith angle range 0° - 45° 
the value of index g is 0.95 for pure proton composition and 
0.91 for pure iron compositon. For energy reconstruction the 
value of g equal to 0.93 was chosen. 



To reconstruct the EAS energy from Cherenkov light flux 
one needs to know an absolute sensitivity of Cherenkov detec- 
tors and the atmosphere transparency. To avoid these problems, 
the method of normalization of the integral energy spectrum to 
the reference one is used. The reference energy spectrum was 
obtained at the QUEST experiment |@]. In that experiment, five 
wide-angle Cherenkov light detectors were installed at Gran- 
Sasso for common operation with the EAS-TOP array. The 
analysis of CORSIKA simulations shows a strict correlation be- 
tween the size/energy (N e /Eo) ratio and the steepness (P) of the 
EAS Cherenkov light lateral distribution. 

The relation between N e /Eo and P (Figj2j is independent 
both from the mass of primary particle and the hadronic inter- 
action model used for the simulation and provides the primary 
energy from the measurement of N e and P. To reconstruct the 
LDF steepness P, the knowlege of PMT absolute sensitivity and 
the atmosphere transparency is not needed. The integral energy 
spectrum of cosmic rays obtained in the QUEST experiment is 
used as the reference one now. The integral energy spectra ob- 
tained for each night of Tunka- 133 operation is normalised to 
that reference one. 

3.3. Reconstruction ofX max 

Recording the pulse waveform for each detector allows to 
use two methods of X max reconstruction, which were developed 
for our experiment. The first is based on the shape of the LDF 
and called P-method. The second method, the W-method, is 
based on an analysis of time width of the Cherenkov pulses. 

The LDF shape is described by an expression with a single 
parameter, the steepness f 071]. Pis strictly connected with the 
distance from the core position to the EAS maximum ]8|]: 

H max =A-BP (2) 

MC simulations shows, that this relation does not depend 
on energy, zenith angle of the showers, mass composition and 
the model of nuclear interaction used for the simulation. To get 
a uniform estimation for P over a wide range of energies, we 
remove from the analysis the detectors at core distances more 
than 200 m during the last step of parameter reconstruction. 



The W-method uses the sensitivity of the pulse width at 
some fixed core distance to the position of the EAS maximum. 
We fixed this distance to 400 m and recalculated the values mea- 
sured with detectors between 200 and 400 m away from the core 
to this distance. To characterize the pulse duration, we use the 
effective pulse with r 



T = S/(1.24 -Ana), 



(3) 



where S is the area under the pulse in [dig.counts * ns], A max in 
[dig.counts] is its maximum amplitude. 

To recalculate the pulse width to 400 m, the width-distance func- 
tion (WDF) is used. This function was constructed on the ba- 
sis of CORSIKA simulation and descibed in J5[]. It was also 
shown 1 5], that the value of t(400) is connected with the thick- 
ness of the atmosphere between the detector and X max (AX max = 
Xq/cosO - X max ) by the expression: 



kXmax - C - D log T 4 00. 



(4) 



This relation is correct for any primary nucleus, any energy and 
zenith angle of the shower and any interaction model, as in the 
case of LDF steepness mentioned above. 

4. Energy spectrum 

4.1. Energy spectrum for inside events 

The data taking by the Tunka-133 array started in Octo- 
ber 2009 and has been continued during two winter seasons of 
2009-2010 and 2010-201 1. As a result, the data were collected 
for 597 hrs of 102 clean moonless nights. Only the nights with 
more than 2 hr of clean weather were taken into account. The 
average trigger rate was about 2 Hz. The number of recorded 
events was about 4 ■ 10 6 . 
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Figure 3: Energy spectrum for two zenith angle intervals, circles - 0° 
30°, squares - 30° < 6 < 45° 
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For construction of the energy spectrum the events with core 
position inside a circle of radius R < 450 m from the center of 
the array were selected. The EAS zenith angle for the events 



collected into the spectrum is limited by 45°. The threshold en- 
ergy of 100% registration efficienty for chosen area and zenith 
angles is 6 • 10 15 eV. The energy spectra obtained for two dif- 
ferent zenith angular ranges (0° - 30° and 30° - 45°) are in 
good agreement for energy more than 6 ■ 10 15 eV (FigfJJ). It 
was found, that the energy spectrum obtained for the first sea- 
son of the array operation (2009-2010) is in good agreement 
with the energy spectrum from the second season (2010-2011) 
(Fig|4]i up to the energy 6 • 10 16 eV. The difference in intensity 
at higher energies seems to be only due to statistical fluctuation. 

The indication of a "bump", found in the first season at 
the energy 8 • 10 16 eV (9J], is not seen for the second season 
spectrum and practically disappears in the combined spectrum 
(FiggJ. 
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Figure 4: Different energy spectra for two seasons of the array operation. Cir- 
cles - 2009-2010, squares - 2010-2011 

The combined energy spectrum for two seasons of the array 
operation is presented in Fig|5] It contains ~40000 events with 
energy E > 10 16 eV and and -400 events with E > 10 17 eV. 
The energy spectrum of Tunka-133 is compared with that of 
Tunka-25 [6], the predecessor of Tunka-133, in this figure. The 
energy spectrum beyond the knee looks rather complicated. One 
can see that the spectrum can be fitted by power laws with 3 dif- 
ferent power law indexes: 

3.21 + (0.01) s ' fl ' + (0.05)*-"' for 6 • 10 15 - 2 ■ 10 16 eV, 
2.95 ± (0.01)"°' ± (0.05) sys ' for 2 • 10 16 - 10 17 eV, 
3.15 ± (0.11)*"' ± (0.05) sys ' for 10 17 - 10 18 eV. 



4.2. Energy spectrum for outside events 

It was found that the energy spectrum including the out- 
side events can be reconstructed too and this spectrum is in 
good agreement with the energy spectrum for inside events be- 
yond some energy threshold. Fig|6]presents a comparison of the 
spectrum for inside events (R < 450 m) with spectra for events 
with core inside circles R < 550 m, R < 650 m and R < 800 m. 

It is seen that the threshold energy of 100% registration 
efficiency increases with increasing of the radius R, but for 
Eo > 6 • 10 16 eV all spectra are in good agreement. Based 
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Figure 5: Summary energy spectrum for two seasons of Tunka-133 operation ( 
red circles) and Tunka-25 energy spectrum ( green circles). 
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Figure 6: Comparision of energy spectrum for inside events (black circles) and 
spectra for outside events R < 550 m ( red circles), R < 650 m (green circles) 
and R < 800 m (violet circles). 



on such the result we construct the combined energy spectrum 
(FigQ from the events with R < 450 m for E < 6 ■ 10 16 eV 
and events with R < 800 m for the higher energy. This com- 
bined spectrum contains about 1200 events with Eq > 10 17 eV. 

The shape of the combined spectrum is nearly the same as 
for inside events spectrum and can be fitted by the power laws 
with 3 different indexes: 

3.21 + (O.OIY"" + (0.05)^' for 6 • 10 15 - 2 • 10 16 eV, 
2.93 + (0.01)'"" + (0.05) !V " for 2 • 10 16 - 10 17 eV, 
3.20 ± (0.06)*"" + (0.05)' 5 "' for 10 17 - 10 18 eV. 



4.3. Comparison with results of other experiments 

A comparison of the Tunka-133 spectrum with the results 
of other experiments is presented in Figj8] Our spectrum is in 
good agreement with spectra of KASCADE(Jlo|]), Tibet (JH) 
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Figure 7: Combined energy spectrum. Green circles - Tunka-25 data. 

and GAMMA ([12]) arrays. The energy range covered by our 
spectrum (10 16 - 10 18 eV) is nearly the same as covered by the 
KASCADE-Grande array data fl3n . 

Both spectra reproduce the same structures: decrease of 
power law index at 2 ■ 10 16 eV and an increase at 10 17 eV. The 
difference in absolute cosmic ray flux intensity for Tunka-133 
and KASCADE-Grande spectra is smaller than 30% even at 
10 17 eV where the difference reaches the maximum value. It 
should be mentioned, that if for g (see expession 1) a value 
of 0.91 is used, the difference in the intensity becomes smaller 
than 10%. 

5. Phenomenological approach 

The parameters derived from CORSIKA simulation may 
slightly differ from the experimentally measured parameters. 
For instance, the linear relation between P and H max observed 
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Figure 8: Tunka-133 all particle energy spectrum in comparison with results of 
other experiments. 
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for MC simulation (equation 2) may hold also for the experi- 
ment, but with a slightly different slope. For our recalculation 
procedure we used the slope derived from experimental data, 
not from MC, i.e choosing a "phenomenological approach". 

In our case we can use the zenith angle dependence of the 
parameter P. This experimental dependence for the fixed en- 
ergy Eq — 10 eV is shown in Fig(9] This mean dependence 
was constructed using all the 14400 events from the energy bin 
16.0 < logio(Eo/eV) < 16.1. The mean zenith angle can be 
recalculated to the mean distance to the EAS maximum. To 
make this recalculation we use the model of the atmosphere 
from 01411 for the real experimental conditions < t >= -30°C 
and Xq = 965g ■ cm -2 . This model gives the following expres- 
sion for the inclined distance to the EAS maximum in km: 



(5) 



where Hq = 96 km. 

To fix the absolute value of < H max > we need to fix the mean 
value < X max > for the above mentioned energy. The most 
reliable experimental estimation of the mean depth of maxi- 
mum by the data of our previous Tunka-25 experiment is < 
X max >= 560g • cmT 2 [8]. Recalculated to the mean < InA >, 
this value is in good agreement with the results of some other 
experiments 11511 . The so derived experimental recalculated de- 
pendence H max vs. P is shown in FigJTOl It can be fitted with a 
linear expression: 



= At - B] ■ P 



(6) 



where Aj = 12.67 km and B\ = 2.09 km. 

This expression was used for the estimation of X max for each 

individual event. 

A similar approach for the second parameter sensitive to the 
position of EAS maximum t(400) can be used but for higher 
energy because of the greater energy threshold of this method. 



The zenith angular dependence of T400 fot the logaritmic energy 
bin 16.5 < logm(Eo/eV) < 16.6 eV is shown in Fig[TT] The 
second method normalized to the X max obtained by the first one 
for the energy 3 ■ 10 I6 eV to fix the absolute scale of AX max . 
This approach results in an expression connecting T400 with the 
thickness of matter between the array and the EAS maximum 
in g ■ cm~ 



-2. 



AX mfl , = Ci - A ■ log t(400) 

where C x = 3493 g ■ cm- 2 , D\ = 1689 g ■ cm- 2 . 
The result is shown in FigfT2l 



(7) 



6. Experiment: < X inax > vs. E 

The experimental dependence of mean < X max > vs. pri- 
mary energy Eq obtained with two methods described above in 
the energy range 5 • 10 15 - 3 • 10 17 eV is presented in Fig[T3] The 
new measurements are compared with that obtained with our 
previous array Tunka-25 and with the theoretical curves simu- 
lated with QGSJET-01 model for primary protons and iron nu- 
clea. The first conclusion is, that the threshold of the W-method 
is higher than that of the LDF steepness method, but the exper- 
imental points obtained by both methods coincide within the 
statistical errors. 

Much higher statistics of Tunka-133 points has led to the 
much smoother behavior of the experimental dependence com- 
pared with the Tunka-25 data. The experimental points go closer 
to the iron curve with energy grow from the knee to about 
10 17 eV. There is a tendency of a backwards movement of the 
experimental points to the proton curve at the energy more than 
10 17 eV, but the statistical errors are too big to insist on such 
conclusion. 

The mean values of < X max > can be recalculated to the 
mean values of < InA > by a simple method of interpolation 
taking into account the corrections to the asymmetry of the X max 
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Figure 9: P vs.9 dependence from experiment. 
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Figure 10: H max vs. P dependence from experiment. 
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Figure 1 1 : T400 vs. 6 dependence from experiment. 



Figure 13: Experimental < X max > vs.Eq dependence for Tunka-133, com- 
pared with Tunka-25. 
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Figure 12: AX maI vs.r(400) dependence from experiment. 

distribution, estimated at our previous work Jj]. The result of 
such approach for the points derived from the LDF steepness 
analysis are shown in Fig[14] We have to note that this proce- 
dure can give different absolute values of < InA > for different 
supposed model of nuclear interaction. The model QGSJET- 
01 we use for the analysis provides the highest position of the 
EAS maximum as compared with the other models used now 
for simulations. The most deep position of EAS maximum can 
be obtained using the QGSJET-II-03 model. The mean differ- 
ence in Xmax between these models is about 20 g ■ cm~ 2 . Using 
of this last model can increase the estimation of < InA > to 
about 0.8 for the same experimental value of < X max >. 

The experimental points are compared with that obtained 
from the analysis of the muon/electron ratio in the KASCADE 
1 1511 experiment in Fig [14] 



7. Other Results 

7.1. Measurements of radio signals from EAS 

To study whether, in addition to the Cherenkov light, also 
the radio emission of air showers can be detected at Tunka, a 
SALLA antenna IU611 was deployed at Tunka cluster 7 in sum- 
mer 2009. SALLA is read out simultaneously with the PMTs 
of that cluster. Indeed numerous radio pulses could be de- 
tected, but most of them are due to RFI by the Tunka elec- 
tronics or other background. However, for Tunka events with 
a high PMT signal, i.e. a high primary energy, an accumulation 
of events with a radio pulse shortly before the PMT pulses is 
observed (see Fig(l5]for an example). Selecting these events, 
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about 70 radio candidate events have been found within the sea- 
son 2009/2010. 

Several tests have been performed to confirm whether the 
detected radio pulses are really from air showers or due to any 
kind of background generated by Tunka. One test was to check 
for which events in a selection of high quality and high energy 
Tunka events a radio pulse could be detected (see FigfToTl. Con- 
sistent with the results of other radio experiments (e.g., LOPES 
|[l7lo . a detection is more likely in the east-west aligned an- 
tenna and for high energies. Furthermore, the ratio between 
the amplitude in the east-west and north-south aligned antenna 
has been compared to RE AS 3 111 811 simulations. For 7 out of 
1 1 events which show a radio pulse in both polarizations, the 
amplitude ratio is compatible within the measurement uncer- 
tainties estimated with a formula derived for LOPES [19]. This 
is another indication, that at least most of the detected radio 
candidates are due to the radio emission of air showers. 

7.2. Time calibration 

For time synchronization of the optical modules of the Tunka- 
133 array, a dedicated calibration system has been developed 
1 20]. The calibration system is based on a powerful nanosend 
LED light source especially designed for the system. The light 
source consists of four high InGaN/GaN royal blue (A max = 450 
nm) LEDs with 100° full emission angle. Thus, the light source 
covers 2n azimuth angle. 

Each LED of the source is driven by its own avalanche 
transistor driver and provides ~ 10 12 photons per pulse, with 
~4ns pulse width (FWHM). All LED drivers of the source are 
triggered simultaneously. Since the light source is installed a 
few meters above the ground it can, in principle, illuminate 
the whole array, given that all optical modules of the array are 
equipped with appropriate reflectors. 
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Figure 16: Energy distribution of a selection of Tunka events and the radio 
amplitude for those of these events which show a radio pulse in at least one 
polarization: radio pulses are preferred visible for higher energies 



8. Plans for upgrading 

The Tunka- 133 facility will be significantly upgraded until 
fall 2013. In 2011 we will install additional clusters of optical 
detectors at large distance from the center of the array ("distant 
clusters"). In 2011-2013, we plan to install radio antennas and 
scintillation detectors for common operation with Tunka-133. 

8.1. Deployment of six distant clusters. 

To improve the accuracy of energy reconstruction for the 
events with the cores located outside the array geomentry we 
plan to install six new clusters at 1 km radius around the center 
of Tunka-133 [9]. These additional 42 optical detectors will 
increase the effective area for Eq > 10 17 eV by a factor of 4. 
The first distant cluster was deployed in autumn 2010, the next 
five will be deployed during summer-autumn 201 1. 
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Figure 15: One of the radio candidate events. A radio pulse is detected in both 
polarizations shortly before the signal of the PMTs. 



8.2. Net of radio anntenas at Tunka 

The positive results of common operation of the first S ALLA 
antenna and the Tunka-133 array suggest that it would be inter- 
esting to install a net of the same type of antennas at Tunka. The 
principle aim of such a radio-array is to show whether radio 
measurements of cosmic ray air showers allow the same pre- 
cision for cosmic ray measurements as Cherenkov light mea- 
surements. In particular, the precision of the energy and mass 
determination of primary cosmic-ray particles will be investi- 
gated. The determination of the precision is a crucial input 
for the design of the next-geneneration large cosmic-ray ob- 
servatories, since radio measurements are not limited to dark, 
moonless nights as other calorimetric detection techniques, like 
Cherenkov or fluorescence light measurements. 

In autumn 2011 we plan to install additional 2-3 antennas, 
further ~20 antennas will be installed in 2012. 

8.3. Scintillation detectors 

The deployment of scintillation counters within the Tunka 
array provides a cross-calibration of different methods of air 
shower parameters reconstruction. 



7 



In 2010, a test variant of a scintillation detector was in- 
stalled. It can operate independently as well as triggered by 
Tunka-133. The mass production of detectors will start at the 
beginning of 2012. We plan to produce 20 scintillation detec- 
tors with 10 m 2 area each. As alternative possibility, it is un- 
der discussion, to move scintillator counters from Karlsruhe to 
Tunka after the shutdown of the KASCADE-Grande array in 
spring 2012. 

Scintillation detectors will be installed with a spacing of 
1 50- 170 m and cover an area of about 1 km 2 . This will permit to 
perform a new QUEST-like experiment [6] and to obtain a point 
of absolute energy calibration at 3 • 10 16 - 10 17 eV in 2012-2013. 
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[18] M. Ludwig, T. Huege, Astrop.Phys. 34 (2010) 438-446. 
[19] EG. Schroder et al. (LOPES Collaboration), NIM A, V. 662, (Supplement 

1) 2012, S238 - S241 
[20] B.K.Lubsandorzhiev et al, Proc of 32-th ICRC, Beijing 2011, ID =1048 
[21] M.Tluczykont et al., Proc. 31th ICRC, Lodz, Poland. arXiv:astro- 
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8.4. HiSCORE station 



The aim of the HiSCORE project 112 111 is to explore the 
gamma-ray skymap beyond 10 TeV with a future large-area 
(10-100 km 2 ) wide-angle (> 0.6 sr) non-imaging cosmic ray 
and gamma-ray air shower detector. 

HiSCORE is a net of detector stations with 4 PMT-channels 
equipped with Winston cones. The station light-collecting area 
(0.5 m 2 ) will be a factor 16 larger than that of an Tunka opti- 
cal detector, resulting in a significantly lower energy threshold. 
The first such station will be installed in fall 201 1 at the Tunka 
site. In 2012-2013, a net of such stations (HiSCORE Engineer- 
ing Array (EA)) will be installed for common operation with 
Tunka-133. The aim of the HiSCORE EA is to prove methods 
of EAS reconstruction before deployment of a full scale array. 



9. Acknowledgements 

This work is supported by Russian Federation Ministery of 
Science and Education (G/C 16.518.11.7051, 14.740.11.0890, 
P681), the Russian Foundation for Basic Research (grants 10- 
02-00222, 11-02-10005, 11-02-0409, 11-02-91332, 11-02-12138) 
Authors acknowledge J.Oertlin (KIT) for his help in the analy- 
sis of the radio-data. 



References 

[1] N.M.Budnev et al. (Tunka Collaboration), Proc. 30th ICRC, Merida, Yu- 

catan,Mexico, 2007, arXiv: 0801.3037, 
[2] N.M.Budnev et al. (Tunka Collaboration), Proc. 31th ICRC, Lodz, 

Poland, 1 (2009) 1069. arXiv: 1003.0089. 
[3] E.E.Korosteleva et al. Proc. 31th ICRC, Lodz, Poland, 1 (2009) 0492. 
[4] N.M. Budnev et al. (Tunka Collaboration), Nuclear Physics B (Proc. 

Supp.), 2009 190: 247-252 
[5] B.V. Antokhonov et al. (Tunka Collaboration), Nuclear Physics B (Proc. 

Supp.), 2011 212-213: 53-58 
[6] E.E. Korosteleva et al., Nuclear Physics B (Proc. Supp.), 2007 165: 74-80 
[7] EAS-TOP Collaboration and E.E. Korosteleva, L.A. Kuzmichev, 

V.V. Prosin, Proc. 28th ICRC, Tsukuba, Japan, 2003 1: 89-91 
[8] N.M. Budnev et al. (Tunka Collaboration), Bull. RAS: Phys., 2007 71 (4): 

474-476 

[9] B.V. Antokhonov et al. (Tunka Collaboration), 25th Texas Symposium 
on Relativistic Astrophysics 2010, published in Proceedings of Science, 
POS (Texas 2010) 138. 
[10] M. Roth et al. (KASCADE collaboration), 28th ICRC, Tsukuba, 1 (2003) 
139 

[11] M. Amenomori et al., Astrophys.J. 678 (2008) 1 165 
[12] A. P. Garyaka et al., Journal of Physics G: Nuclear and Particle Physics, 
35,115201 (2008) 



8 



